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Abstract: In this paper we address two questions: how energy conservation affects the solution to the non- 
linear equation, and how impact parameter dependence influences the inclusive production. Answering the 
first question we solve the modified BK equation which takes into account energy conservation. In spite of 

t^- the fact that we used the simplified kernel, we believe that the main result of the paper: the small (< 40%) 
suppression of the inclusive production due to energy conservation, reflects a general feature. This result 
Jeads us to believe that the small value of the nuclear modification factor is of a non-perturbative nature. In 

O the solution a new scale appears Qf r = Qs exp(— l/(2as)) and the production of dipoles with the size larger 
than 2/Qf r is suppressed. Therefore, we can expect that the typical temperature for hadron production 

^ is about Qf r (T Q/ r ). The simplified equation allows us to obtain a solution to Balitsky-Kovchegov 

^ equation taking into account the impact parameter dependence. We show that the impact parameter (6) 
dependence can be absorbed into the non-perturbative b dependence of the saturation scale. The solution 
of the BK equation, as well as of the modified BK equation without b dependence, is only accurate up to 
±25%. 
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1. Introduction 



High density QCD that has a long history [1-6], it is actually based on the idea that before summing 
the next-to- leading order correction to the BFKL Pomeron [7], we need to deal with different kinds of 
corrections. Indeed, the scattering amplitude in the leading order BFKL approach is proportional to 

iV BFKL oc o? s g*™*™ (1.1) 

where A B ^p KL oc as ■ We can see two kind of corrections to the amplitude of Eq. (1.1): the first one is the 



next-to-leading order correction to the BFKL kernel that gives Abfkl = otsC\ + a s C2 where C\ and C2 
are constants; and the second corrections due to the exchange of two BFKL Pomerons, namely, 

^2 BFKL Pomeron exchange a ^BFKL K a s s 2Abfkl (1.2) 
Comparing Eq. ( |L~l| ) and Eq. ( |1.2[ ) one sees that for a wide range of energies 

y = Ins > — In ( — J (1.3) 

as \OisJ 

One has also to include the exchange of two BFKL Pomerons that leads to non-linear corrections in the 
evolution equations. 

For these energies the NLO corrections to the BFKL kernel are small . They becomes essential for 
higher energies for which a s y ~ 1 or y = Ins > 1/as- Therefore, for energies in the kinematic region 

V4>y>;|m(-L) (i.4) 

we have to include multi-Pomeron exchange and can neglect NLO corrections to the BFKL kernel. 

Explicit calculations of the NLO corrections to the BFKL kernel [8], show that a s corrections to 
Abfkl are rather large, and have to be taken into account for any realistic and practical estimates. In 
terms of Eq. ( |1.4| ) the large NLO corrections mean that the kinematic range for the high density QCD 
is narrow. Fortunately, we know the NLO corrections to the linear BFKL equations [12] and to the non 
linear (Balitsky-Kovchegov) equation (see Ref. [13] for review and references) quite well. We also know 
that such corrections change the value and energy dependence of the saturation scale [14,15]. Therefore, 
it is necessary to include the NLO correction to obtain reliable estimates of the non-linear effects. 

It is not sufficient to solve the Balitsky-Kovchegov (BK) equation because together with the NLO 
corrections to the Balitsky-Kovchegov equation, we need to include in the calculating procedure the cor- 
rections referred to as Pomeron loops [10,16-22]. In spite of the fact, that we have learnt a lot about these 
corrections, no closed equation exists and the whole procedure turns out to be so complicated, that there 
is no hope of obtaining reliable estimates in the foreseeable future. 

In such a situation the only thing that we can do, is to try to study the influence of different terms in 
NLO corrections to the Balitsky-Kovchegov equation, with the goal of finding the essential terms. After 
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that we can try to estimate the influence of the Pomeron loops and find out what is more important: loops 
or NLO corrections. We are presently at the first stage of such a study. 

The clear example of such approach is the effect of running QCD coupling constant. The running QCD 
coupling constant is believed to be one of the most characteristic features of QCD, and it is important 
to find out how this phenomenon affects the value of the scattering amplitude. Fortunately, the form of 
the BK equation is known [23] and the numerical solution [24] shows that the solution with running QCD 
coupling constant is quite different from the one with fixed coupling. 

The second example is the solution to more general JIMWLK equation which includes all 1/N C correc- 
tion to Balitsky-Kovchegov equation (see Ref. [25]). It turns out that the difference, between the solution 
of the JIMWLK equation and the BK one, is extremely small. 

In this paper we wish to investigate the influence of the conservation of energy on the solution to the 
BK equation. As has been discussed the important observable for the interpretation of the nucleus-nucleus 
collision is the energy loss [28]. In LO the non linear equation is written neglecting the energy loss. There- 
fore, we need to generalize the BK equation in such way that it would respect energy conservation. Such a 
modified equation will allows us to discuss the value and energy and transverse momentum dependence of 
the nucleus modification factor (NMF). Since for single inclusive production the kt factorization has been 
proven (see Ref. [29]) the value of the NMF, that is originated from short distances, is determined by the 
BK equation. The generalized BK equation has been suggested in Ref. [30] and for the completeness of 
presentation we will discuss this equation in the next section. The main section of the paper is the third one 
in which we present the solution to this equation closely following the method proposed in Ref. [31]. We 
find no considerable effects due to energy conservation, and show that the dipole scattering amplitude in 
the region r > 1/Q S leads to the inclusive gluon production that is 20 - 40% smaller than that is predicted 
from the LO BK equation. 

The second problem that we discuss in this paper is the impact parameter dependence of the solution 
both to Balitsky-Kovchegov equation and to the modified version of this equation. It is well known that b 
dependence is one of the most difficult and challenging problems in QCD [32]. The massless gluon leads to 
power-like decrease of the scattering amplitude for large values of b which results in the power-like increase 
of the interaction radius [32] . Such behavior of the interaction radius indicates that the main contribution 
to the solution of the non-linear equation stems from large values of b making the entire approach, based 
on perturbative QCD, inconsistent. The attempts ( only two [33] as far as we know) to solve BK equation 
numerically taking into account b dependence, confirm these pessimistic expectations. The simplified 
version of BK equation allows one to obtain the analytical solution including 6-dependence. It turns out 
that the non-perturbative b dependence can be absorbed in the saturation scale. This observation gives 
us a guide for the solution of the impact parameter problem in high density QCD, but how to include 
this observation into general BK equation still remains unsolved. In section 4 we discuss the inclusive 
production for the solution with impact parameter dependence. 
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2. The modified non-linear equation 

The non-linear (BK) equation has the following form for the dipole scattering amplitude N (r, Y; b) 

9N ^' b) = ^ J d 2 r , K[ry) f 2N h Y ;b- ^(r-r')) - N (r,Y;b) 

- N (r',Y;b-^(r-r')\ N (r-r',Y;b-^r')^ (2.1) 

where Y = ln(l/x),as = FCts , x denotes the Bjorken variable for the dipole-target scattering, and b is 
the impact parameter for the reaction. Kernel K (r, r') in the LO is equal to 



K ( r > r ') = ( f _ r -v )2r ,2 ( 2 - 2 ) 
2.1 The linear equation in NLO 

The linear part of Eq. ( |2.lD can be easily written in the double Mellin transform 

N{d = Hr 2 /R 2 ),Y = Hl/x);b) = J ^^L N (u;,r,b)e" Y + ^ (2.3) 

It has the form 

LoN(u,r,b) = (asXLoil) + a 2 s XNLo{l)) N{u),r,b) or u = lo(j) = a s XLo{l) + o? s xnlo{i)- 

(2.4) 

F° r Xlo(i) we have the well known expression [7]: 

Xlo(i) = 2^(1) - V(7) - V>(l-7) = - + -T— + Xlo(i), (2-5) 

7 1 - 7 

with ip = dlnr(7)/c?7 and V is the Euler Gamma function. We will discuss the form of Xnlo(i) later. 
Xlo(^) = 2 "0(1) — V'(l + 7) — V'(2 — 7) denotes the contribution of the higher twist which we will discuss 
below. 

The NLO BFKL kernel is known [12,26], following Ref. [30] we use the simplified form of the NLO 
kernel: 

• In Ref. [15] the following form of the NLO BFKL kernel is proposed 

asXNLoh) = + - z - ujxlo 7), 2.6) 

V 7 7 1 - 7 + oj 1 — 7 / 

where, 



A(u>) = -11/12 + 0(cu) + n F (^jL- P Gq (u;) P qG {u>) - 1/3) , (2.7) 
with P(oj) being the DGLAP kernel. 
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The singularities in Eq. ( [2.6D describe the different branches of evolution corresponding to the sizes 
of the interacting dipole (or the transverse momenta of partons). The pole at 7 = corresponds 
to the normal twist- 2 DGLAP contribution, with the ordering in the transverse parton momenta 
Q > . . . kij > h,i+i > • • • > Qo, where Qo is the typical virtuality of the target Qo « l/R. The pole 
at 7 = 1 in Eq. ( |2.6| ) corresponds to inverse kt ordering (k^t < ktj+i < . . . Qo)- The other poles, at 
7 = — 1, —2, ... (7 = 2, 3, ... ), are the higher twists contributions due to the gluon reggeization. 



• As can be seen in Eq. (2.6), the main changes in the NLO kernel to the ordinary DGLAP evolution 
are introduced, so as to account for the DGLAP anomalous dimension, and to the part of the kernel 
that describes the inverse evolution. This branch of evolution is moderated by the non-linear effect, 
rather than by the NLO corrections. Indeed, it was shown in Ref. [31] that the term 1/(1 — 7) leads 
to exponentially small corrections in the saturation region. Consequently we can neglect changes in 
the inverse evolution, and keep the BFKL kernel without a shift from 7 to 7 — to. 



• We use the observation of Ref. [34], according to which Eq. (2/7) for A\{oj) can be approximated 
to an accuracy (> 95%), by A\{uj) = 1. 

Hence, we can rewrite the NLO BFKL kernel in a very simple form: 

XNLoil) = -UXLOU)- ( 2 - 8 ) 

Using Eq. ( |2.8[ ) we obtain the full kernel for the linear equation in the form: 

w ( 7 ) = asxioil) (1 - w) (2.9) 

This kernel imposes energy conservation (see Ref. [34]), and describes the NLO BFKL kernel. It does 
not include the contribution coming from inverse ordering, which should be suppressed in the solution to 
the non-linear equation. 

Note that substituting \lo = olsH yields 

7H = a s (~ ~ l) , (2.10) 

which approximates the DGLAP anomalous dimension in leading order so well, that the difference between 
this anomalous dimension and Eq. ( |2.10| ) turns out to be less than 5% [34]. 

Eq. ( 2,10j ) corresponds to Eq. ( |2.6| ) where A\{to) = 1, and the high twist terms are neglected. 



2.2 The non linear equation in NLO 

The main idea of Ref. [30] is to use Eq. (p^) as the kernel for the non-linear equation. 



First one can see that Eq. (23) can be rewritten in the coordinate space as 
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dN (r, Y; b) C F a s 



dY it 



2 



J d^K LO ( r y){i - A} 



(^2N (r',Y;b - ^(f-f')\ - N(r,Y;tj\ (2.11) 

As is shown in Refs. [3,21,27] the non linear evolution equation has a very simple probabilistic inter- 
pretation and can be written as the equation for the Markov chain or, in other words, can be viewed as 
the process of the birth (non linear term) and death (linear term) of the colorless dipoles. Assuming that 
this interpretation is correct in NLO order, we obtain the following equation [30]: 

Ml ^/^ lo(r y ){l _^} „ 12) 

(^2N (r',Y;b - ^(r-r')\ - N (r,Y;bj - N (r',Y;b- ^ (f-r') \ N fr- r',Y;b- \r'X\ 

The advantage of this equation is that energy is conserved by the non-linear term, as well as by the 
linear one. The linear term has been discussed above. We need only to rewrite the non-linear term in 
Mellin transform ( see Eq. (|2.3| )). namely, 



N (r',Y-,b-^(f-f')\ N (r-f',Y;b-^f'\ = (2.13) 

to understand why this term conserves energy, as the derivative with respect to Y leads to extra factor uj 
in the non-linear term. 

3. Solution to the modified B-K equation 

3.1 Simplified kernel 
3.1.1 General approach 

As was suggested in Ref. [31] we replace the full BFKL kernel in LO by the following simplified kernel 

i for z = In (r 2 Q 2 s ) < ; 



ulo (7) = «s < 



7 

(3.1) 

for z = ln(r 2 Q 2 s ) > 0; 



This kernel leads to the leading twist contribution and gives the natural generalization of the DGLAP 
equation that includes two different kind of logs: for z < it sums (agin (r 2 AQCD)) n while for z > the 
new type of appears, namely, [as In (r 2 Q 2 )) n . 
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For the NLO kernel we assume that for z = In (r 2 Q 2 ) < we have Eq. ( |2.10| ) which can be rewritten 

as 

UNLO (7) = otS — -rr- (3.2) 

7 + a S 



This form of the kernel shows that in the NLO, we sum, in addition to [as In (1/ (r 2 Aqcd)) J terms, also 
non-logarithmic corrections. Recall that log contribution stems from the singularity (I/7) in the kernel. 

For z = In (r 2 Q 2 ^ > we parameterize the NLO kernel as the solution to the following equation: 



unlo(i) = a s (1 - ^nlo{i)) { 1 - - + as «} (3.3) 



where we determine the constant k from the condition 



UNLO ( 7 7>>r > 7cr ) = ^NLO (7 7<7 " > 7cr ) (3.4) 



Therefore, to find k it is necessary to discuss the critical value of the anomalous dimension. 



3.1.2 Critical anomalous dimension and saturation momentum 



Using Eq. ( p.l[ ) one can calculate the critical anomalous dimension. It is well known that for this calculation 
one only needs to know the kernel for z < 0( see Refs. [1,14,35,36]). The equation for the critical anomalous 
dimension (7 cr ) has the following form 

duj(-fcr) = Cj(7cr) , g g v 

O^fcr 1 - lex 



Inserting Eq. ( p.l|) in Eq. (3J5) one obtains 



lev = 2 (1 - a s ) (3.6) 
he equation for the saturation momentum has the form [1,14,35,36]: 

HQVQl) = 77^72^ (3-7) 
(1 + as) 

where Y = ln(l/x), and Qq is the scale associated with the initial condition at low energies. 

In the vicinity of the saturation scale z — > 1 the behavior of the dipole amplitude has the form [14,37] 

N(Y;r) oc {^Q 2 ) 1 '^ = e {1 ~^ )z (3.8) 
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3.1.3 The NLO kernel for 7 > j c 

Using Eq. (j3lj) we can calculate 



/ 7<7cr \ 2« S 

Watlo 7 > Ict = , , _ (3.9) 

/ 1 + as 



and Eq. (|3.4|) leads to the following value of k: 



l-a 2 s 



4«5 

fc — - 
1 

Finally, the NLO kernel has the form: 

W7VLO (7) = OS < 



At = - ; 2 (3.10) 



^- for z = ln(r 2 Q 2 ) < ; 



7 + a s 



1 7 + a a « , 2 2 \ > q 



(3.11) 



It is worthwhile mentioning that the kernel of Eq. ( |3.11| ) describes the LO anomalous dimension of the 
DGLAP equations for z > 0, and it takes into account energy conservation. Therefore, we consider it as a 
model that can lead to a reliable predictions. 



3.2 The simplified equation in the coordinate space 
3.2.1 z < 

In this kinematic region we can simplify K^q (r, r') in Eq. ( |2.12| ) in the following way [31], since r' 3> r 
and |r — r'l > r 

j d 2 r'K LO (r,r') 7rr 2 ^ (3.12) 

Introducing n (r, Y; b) = N (r, Y; b) jr 2 we obtain 

d 2 n(r,Y;b) r 3 



drMn(l/(r 2 A 2 )) 1 ^ 



«s{l - ^7} (2n(r,F;6) - r 2 A^ D n 2 (r, Y; b) ) (3.13) 



3.2.2 z > 

The main contribution in this kinematic region originates from the decay of the large size dipole into one 
small size dipole and one large size dipole. However, the size of the small dipole is still larger than 1/Q S . 
This observation can be translated in the following form of the kernel 

d 2 r'K LO (r,r') ^ vr / % + vr / ( 3 - 14 ) 

Ji/Ql(Y,b) r' 2 A/Q2(K,6) {r-n 2 
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One can see that this kernel leads to the In (r 2 (5;?)-contribution. Introducing a new function N (r, Y; b) = 

2 

J r dr 2 N (r,Y; b) jr 2 one obtain the following equation 

d 2 N(r,Y;b) r d X (( 8 dN(r,Y;b) \ - \ 

This equation does not depend on the value of the saturation momentum. Because of this we can assume 
that the solution of this equation depends on one variable z = In (r 2 Q 2 ^j . The result that the solution in the 
saturation region depends on only one variable, has been proven (see Ref. [35]). The value of the saturation 
momentum depends on two variables: on rapidity Y and impact parameter b. Since the equation does 
not depend explicitly on Q s and b, the only information that we have on these variables stems from the 
matching with the equation at z < on the critical line. Transforming from variable Y and r in Eq. ( |3.15 ) 
to variable z = In {r 2 Q 2 s ) we have 

d 2 N(z) _ {l + asf ( x 4a s d dN(z)\~ T/ ^ ^ d f dN {z) 



dz 2 4 I 1 — a| dz dz 



It worthwhile mentioning that the BK equation in this kinematic region has the form 
3.3 Solution to BK equation with the simplified kernel 



For the sake of completeness we start by recalling the solution to BK equation in the form of Eq. ( |3.17 ), 
obtained in Ref. [31]. Introducing 



N(z) = J dz' (l - e -^')) (3.18) 
Eq. ( |3.17| ) can be rewritten in the form 

»£> = i f « (1 - e-*'») ( 3, 9) 
Considering dcfi(z') /dz = D (0) we can rewrite Eq. ( |3.19| ) as 

Differentiating Eq. ( ^j.20|) with respect to <j> we have 

1 - e-f) (3.21) 



dD{<j>) _ 1 1 
dd> ~ 4 Die 
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which leads to 



and finally, 



D 1 



+ e" 



V2 



'4> V<A + e"* - 1 
The boundary conditions at z = have the form [1,31,37] 

0(z) ► (/> < 0(z); 

din (0 (z)) z ^q-^ i z-K)+ din (0 (z)) _ 
dz 2 dz 

Eq. ( |3.25| ) follows from Eq. ( |3.8[) for the kernel given by Eq. (3.1) 



(3.22) 
(3.23) 



(3.24) 
(3.25) 



3.4 Solution to the modified BK equation 

Using Eq. ( gig ) we can reduce Eq. fl3,16| ) to the form 



e"*W 2V (z) + " S (1 ^ - a 5 (V^ N (z) 



or 



/ (1 - a s N 



1 - a s 



N(z) + a s (1 - e 



-0(2) 



(3.26) 



(3.27) 



Considering being function on N and neglecting the terms that are proportional to a|, we can rewrite 
Eq. (ET27D as 



1 



(l + «5) 5 



AT 



1 - a 5 TV v ^ 1 - exp <| - ( iV 
The BK equation for cj) [N] has the form 



1 



} 



N 



+ a s (l- exp{-</>(iV) })) (3.28) 



JV 



(3.29) 



4 1 - exp { - (\> (n\ } 

In comparison with the BK equation Eq. ( 3.28; ) (see Eq. ( |3.29| )) this contains a new feature: 4>'^ — > oo 



at N — > 1/as- In the vicinity of N = 1/ag we can simplify Eq. (|3.28 ), namely, 

1 /{l + a s ) 2 N y 



N 



1 - a s N 



I - exp{-0(iv)} 



(3.30) 



since the first term in the bracket in Eq. ( 0.28 ) is of the order of > 1/as- 
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The solution to Eq. ( 3.3C| ) has the form 



+ e" 



(l + a s ) 2 1 



In ( 1 — asN ) + as 



(3.31) 



From Eq. ( p.31| ) we can see that (ft is large in the region where asN <C 1. Having this in mind we can 
replace Eq. (|3.28|) by the following equation 



1 



1 - a s N 



N + a s 



(3.32) 



which gives 



N) = I -a s 



(l + «s) 5 



N 



(3.33) 




Figure 1: Dependence of 



on N 



for solution to Eq. ( fHgj ) (<f>W and to 
Eq. O (0(°). a s = 0.2>o = 0.3, 
N = 0.6 



ion 



In Fig. we illustrate the behavior of (ft (^N^ for the solut 
to Eq. ( 3.28| ) and Eq. ( p.29| ). The difference between modified 
equation and BK one is clearly seen, and it is well reproduced by 
the simplified solutions of Eq. ( 3.31| ) and Eq. ( [3.33 ). However, 
for small values of N we have to solve equation numerically to 
determine the accuracy of the estimates. Having (ft we can 
find N as function of z. Indeed. 



1 — exp ( —(ft ( N 



dN(z) 
dz 

and the following equation gives N as function of z 

dN' 



r-N 



No 1 - exp [-(ft [N 



(3.34) 



(3.35) 



The value of N (z = 0) = A^ can estimated from Eq. ( |3.8|) . Indeed, 
in Ref. [31,37] was shown that for z < but in the vicinity of z = N = f*^ dz' (l— exp (—(fto exp ((1 — j C r)z)) 
. Taking this integral we obtain the value for iVo- 

In Fig. | the solut ions for modified equation and for the BK equation are plotted. The comparison 
between them is illustrated in Fig. || in which the ratio N^/N^ is shown, where iV"W and N^ ' are 
solutions to Eq. ( |3.28|) and Eq. ( 3.2S| ), respectively. One can see that the effect could be as large as 20% 
(see Fig. |) 
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Figure 2: Dependence of function <j>(z) (Fig. |^-a) and the scattering amplitude N(z)( Fig. ||-b) versus z for solution 
to Eq. ( Uf ) (AfW) and to Eq. (JV^). a 5 = O.2,0 O = 0.3, A> = 0.6 



3.5 Inclusive production 

An interesting case is to compare the inclusive production for both equations. The inclusive production can 
be calculated using kt factorization, this was proven for the case of hadron-nucleus interaction in Ref. [29]. 



da 



dy d 2 pT 



(3.36) 



2C F 1 



a s (2vr) 4 p 2 T 



J d 2 bd 2 Bd 2 r T j^ V 2 T N G {y l = \n{l/x l );r T -b) V 2 T N G (y 2 = ln(l/x 2 ); r T ; \b - B|) . 



In the saturation region S/ 2 Ng = V 2 (2N — N 2 ) oc exp (—2cj)(z)) and, therefore, 

da 



dy d 2 pT 



oc e 



(3.37) 



Using Fig. ^-a we can see that the effect is large, and practically for z > 3 the ratio of exp (— 4</>) is negligibly 
small. In other words, the production of the gluons with transverse momenta smaller than Q s exp[— 1.5) is 
small. 

To estimate the inclusive production we need to know V 2 N G - It is easy to see that V 2 iV G = 
(l/r 2 )d 2 N G /dz 2 . Eq. (|]27|), Eq. (&2l\ ), Eq. (|]|§ and Eq. (|3T35| ) allow us to calculate <j)' z and 4>" zz . 
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Indeed, 

d#(°)(z) 



dz y/2 



_L ^(o)( z ) + e -0(°)( 2 ) _ i; (3.38) 



dz 2 4 V ' )' 



dz 1 - a s N(z) 

- a 5 iV(z) N2 



(3.39) 



(3.40) 



Q±*s£(i _ ^ (1 »W) _ 2d g ^ W(z) e ~* (1) (') fl- e -* (l) W 



+ ^ ^ = ^ '-\ 3.41 

1 - a S N(z) 




Figure 3: Dependence of the ratio N^/N^ versus z (as = O.2,0 O = 0.15, N = 0.6 

In Fig. [2] the solutions for modified equation and for the BK one are plotted. The comparison between 
them is more illustrative. 
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For y\ = y<2 Eq. ( 3.36j ) can be re-written in the form 
da 2C F 1 



dycPpT as2(2-K) 



d 2 b d 2 B 



+00 



dze z Jo[e2 z x± 



d 2 N G (z;b) d 2 N G (z;\b-B 



dz 2 



dz 2 



. (3.42) 



The integral over z in Eq. ( gjg ) includes the region for z < 0, where the function </>(z) behaves as is 
shown in Eq. (|3,8|), however the kinematic region where we can trust Eq. fl3.8| ) is rather narrow [37]. Using 
Eq. ( |3.1| ) and Eq. ( |3.11| ) we can find the full solution for (j){z) at z < 0: 



(z;z < 0) = (/> exp (L + z) - L - z) 



;z < 0) = <p exp ((1 + a s ) V L ( L + z ) ~ (1 + «s) + 



(3.43) 
(3.44) 



where xj_ = pt/Qs (Y) and L = 4a5y for BK equation and L = 4ag/(l + as) 2 y for the modified BK 
equation. 

Using these expressions we can calculate d 2 N/dz 2 for both cases at z — > for z < 0: 



d 2 N(°) z ^ - fo( ^0+ d 2 N^ 

dz 2 y 4 I + Lj * 4 ^ dz 2 



dz 2 



(1 + as) 



1 + 



;i + as) z^o+ d 2 N^ 
4 < dz 2 



(3.45) 
(3.46) 



From Eq. (|3.45|) and Eq. ( ^46|) we have that the double derivatives at z = only match at large values of 
L. The calculations for the ratio 

da (i) d(J (0) 

dy d 2 px ^ dy d 2 pT ^ 
are shown in Fig. [|. One can see that the energy conservation could lead to the suppression of the inclusive 
production by a factor of two. 



3.6 New scale for freeze-out 



The most striking result of the modified BK equation is the appearance a new scale. Indeed, Eq. ( |3.32 ) 



and Eq. (3.33) show that for z when asN — > 1 , the value of increases. Formally speaking for asN > 1 
<fi becomes an oscillating function which gives negligible contribution to the inclusive cross section (see 
Fig. H) . Since in vicinity asN ~ 1 the value of <j> oc l/a s > 1, we can conclude from Eq. ( 3.18 ), that 
N = z. Therefore, the equation asN = 1 leads to z = l/as, which gives a new scale 



Qfr 



Qs exp 



2a s 



(3.48) 



The physical meaning is clear: the dipoles with the size larger than 2/Qf r will not produced. In other 
words, in the collisions we produce dipoles with the size 2/Qj r with the typical temperature T = Qt r . 
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Figure 4: Dependence of d 2 N/dz 2 versus z for the solution to BK equation (dashed line) and to modified BK 
equation at different values of as. 

It should be stressed that we reconstruct the form of the non-linear term in the modified BK equation 
assuming the probability interpretation of the parton cascade, which was proven for the LO BK, however 
we believe on physical grounds that it will also be correct in the NLO. However, the NLO corrections also 
lead to the transition one dipole decays into three dipoles and induced by this transition term in the linear 
part of the equation. We neglected these terms since they induce the A^ 3 term in the equation. This term 
is small in the vicinity of the critical line ( N oc a$ on it [31]) and only this kinematic region contributes 
to inclusive production. 

The failure to obtain large corrections from energy conservation forces us to believe that non-perturbative 
contribution will be the reason for the small value of the NMF( nuclear modification factor). It should 



be stressed that the dipole rescatterings has been taken into account in Eq. ( 3.36 ). A separate question 



is, could we trust the kt factorization formula of Eq. ( p.36|) in the case of nucleus-nucleus collisions? The 
rigorous answer is no since there is no proof of the kt factorization for this case. However, it is difficult 
to believe that the violation of the kt factorization will be responsible for the value of NMF, Indeed, in 
Eq. ( p6| ) and Eq. flOg ) the main contribution stems from the vicinity of the saturation scale (the largest 



one in nucleus-nucleus scattering) and in the case of the two scale problem: high pr and the saturation 
scale, kt factorization has been proven in Ref. [29]. 

We think that the most interesting outcome of this solution is the appearance of the new scale 
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Qfr = Qs exp(— 1/(20:5)) which characterize the temperature of the produced hadron ( T « Qf r ) . 
The production of dipoles with the size larger than 2/Qf r are suppressed. Much more work is need to 
answer the question whether this scale is an artifact of our simplification or a general feature of the energy 
conservation. 



4. Impact parameter dependence 



In this section we address the problem of the impact parameter (6) dependence for both BK equation and 
modified BK equation. It is well known that b dependence is one of the most challenging questions in 
perturbative QCD [32]. Indeed, in QCD we have the massless gluon, and because of this at large values 
of b we are doomed to have power -like decrease of the scattering amplitude in both linear and non- 
linear equation. Such a decrease results in power-like increase of the radius of interaction R oc s A , which 
indicates that the large values of b are dominant in the equation, which has been derived in perturbative 
QCD [32]. Therefore, at first sight the approach based on non-linear equation appears to be inconsistent. 
Unfortunately, as far as we know there exist only two papers (see Ref. [33]) where the numerical solutions 
were found to BK equation including b dependence. These solutions confirm the result that large values of 
b become essential and some modifications of the approach based on the BK equation, are required. 

In such a situation the analytical solution that includes b dependence even for the BK equation with 
simplified kernel, can be instructive. 



4.1 Solution with the impact parameter dependence 

For simplified kernels of Eq. ( |3.1| ) and Eq. ( |3.11| ) the 6-dependence can be easily found following Ref. [31]. 
Indeed, in both equations the kernel for r 2 Ql < 1 coincides with the DGLAP kernel. In DGLAP evolution 
we sum the log of transverse momenta of the partons. In this integral the momentum transfer q enters only 
as the low limit of integration over transverse momentum (say k±). The integral has the following form 

HQ 2 /* 2 ) = ( 4 -D 



However, the low limit in Eq. (4.1) is q only in the region where q > Qq where Qq is the scale from with 
which we can apply perturbative QCD. Therefore, in the kinematic region where q < Qq no dependence of 
q appears in the QCD evolution. However, the non-perturbative corrections generate the non-perturbative 
form factor, and the general form of the solution for say, gluon density, has the following form 

N (z < 0; r, b; Y) = 5(b) N (z < 0;r;Y) ^> S (b) e^ 1 " 7 ^ z (4.2) 

In the vicinity of the saturation scale Eq. ( |4.2D together with Eq. ( |3.8[) leads to 

N(z^0;r,b;Y) = (r 2 Q^Y))" S (b) = (r 2 Q^Y-b)) 1 '^ (43) 
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x ± = p x /Q s 

Figure 5: Dependence of the ratio R = J^ip T I dy'd i p T ^ or V 1 = V 2 versus = Pt/Qs and different values of L 
where L = AasY for the solution of BK equation and L = 4as/(l + as) 2 Y for the solution of the modified BK 
equation. In the picture the value of L for the BK equation is shown, (as = 0.2, <fio = 0.15, No = 0.6) 



where 

Q 2 s (Y;b) = QliY) S^(b) (4.4) 

The kernel of the BK equation as well as the modified BK equation does not depend on b. Therefore, 
the only information about 6-dependence stems from the 6-dependence on the critical line which can be 
absorbed in the dependence of the saturation scale on b as is given by Eq. (fO|). Finally, the solution of 
the equation will be the same as without 6-dependence but with new variable 

z(z;b) = In (r 2 Q 2 s (Y;b)) = z + ]n(S (b)) (4.5) 

V J 1 — 7cr 

This new variable alters the estimates for the new scale for freeze-out as one can see from Fig. g that 
shows N as function of z at different values of b. One can see that the value of z at which N — > 1/a.s is 
different for different b. In terms of ion-ion collisions this dependence can be translated to dependence on 
the centrality cuts leading to temperature dependence for different centrality. 



4.2 Inclusive production 

Eq. ( p.42| ) can be re-written in the form (for y\ = y^] 

da 2C F 1 [ + °° , _ z T / i 2 \ r f d 2 N G iz(z;b) 



dycPpr as2(2ir) 3 x 



3 rJl ' 



J dze~ z J (e2 z Xl _) [y d 2 b Lj (4.6) 
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Figure 6: N versus z and for solution to modified BK equation at different values of b ( see Eq. ( |4.7D ). The dotted 
line shows N = 1/cts for as — 0.2 

where x± is defined in Eq. fl3.42| ). 

As we have mentioned in section 3 we calculated the inclusive cross section assuming that N oc 
(R — b) Here, we choose the exponential dependence for S (b) 

S(b) = exp(-6 2 /^ 2 ) (4-7) 

We take R 2 = UGeV~ 2 for the nucleon and for nuclei we choose R 2 = (2/h)R 2 vs where Rws is the 
radius of the nucleus in Wood-Saxon parameterizations (Rws = 6 A 1 / 3 GeV^ 1 ). 

Changing variable b —> p = bj R oneobtains 



da 2C F 1 



d 2 N G (z(z;p) 

dze- z J [e^x ± ) \ I dp 2 L \ = T ( x ±) ( 4 - 8 ) 



- i 



(ttR 2 ) 2 dyd 2 p T a s 2(2vr)3 x^V-oo V / \ J dz 2 



Function T (x±) does not depend on the target. Therefore, Eq. ( |4.8| ) shows that for the scattering of two 
identical nuclei at y\ = yi we expect the scaling behavior: the inclusive production depends only on variable 
x±, instead of dependence on the number of nucleons (^4), rapidity y = y\ = yi and pt- 

Fig. [?] shows the influence of b dependence of the inclusive production. 

In Fig. the ratio 



R(b) = T(with b dependence) / T[ without b dependence) (4-9) 

is plotted. 
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Figure 7: The ratio R (b) ( see Eq. ( |4.9| )) versus x± for the solution of the modified BK equation (curve 1) and for 
the BK equation (curve 2). Curve 3 gives the ratio of T{x±) for BK equation to the modified BK equation without 
6-dependence while curve 4 gives the same ratio taking into account 6-dependence. 



One can see that that the solution of the BK equation with b dependence, differs from the solution 
of the same equation without taking into account the 6-dependence, by + 25% at small values of x± and 
and by -25% at large values of x±. It is surprizing that the solution of the modified BK equation leads to 
the same inclusive production at large xi both for the cases with b and without b dependence. However, 
at small x± the solution with b dependence exceed the one without 6-dependence by 35%. The differences 
between solutions to modified BK and BK equations are shown in curves 3 and 4 in Fig. [|. The solutions 
with 6-dependence could differ by a factor of three at large values of x± . 

T (x±) as a function of x± is plotted in Fig. ||. Note that the nuclear modification factor(NMF) can 
be easily re-written through T(x±), namely, 

d 2 N AA 

NMF = — = -4- j ^ j . (4.10) 



N ^ Sfe " A2/3 T[x ± 



!a,A 



where N co u is the number of collisions. Using the calculated T (x±) (see Fig. |8|) and the fact that 
Qs,a/Qs,n = A 1 / 6 we see that for the gold NMF = (0.2 0.3) x 1.7 =0.3 0.5 for x± = 0.1 -=- 10 ( 
pj_ = 0.14-^14 GeV) which is in about 2 times larger than the experimental NMF ( see RHIC data [38-41]). 
However, the ratio of the radii (R\/ 'R% = 1.7) as well as Gaussian parametrization cannot be considered as 



- 19 - 



reliable. Therefore, we can conclude that both BK equation and the modified BK equation lead to rather 
small NMF but perhaps, this suppression is not sufficient to describe the experimental data. 




Figure 8: T(x±) (see Eq. fl4.8| ) versus xj_ = pj_/Q s for the solution of BK equation that takes into account 
6-dependence. 



5. Conclusions 



In this paper we study two questions: the influence of energy conservation on the solution to the non-linear 
equation; and the 6-dependence of the solution. The influence of energy conservation on the solution of 
the BK equation was investigated taking the impact parameter dependence in the form of (b — R). The 
answer to this question is that the energy conservation reduces the value of the inclusive cross section by 
20-40% at reasonable values of the QCD coupling (as ~ 0.2). These estimates are based on the modified 
BK equation (see Eq. ( 2.11[ )) and on the form of the kernel given in Eq. ( 3.1 1| ) , which were derived from 
the form of the NLO corrections to the BFKL kernel, given in Eq. (2^). The derivation of Eq. (|2.8|) stems 



from two observations: (i) the fact that Eq. ( |2.10| ) describes the exact value of the anomalous dimension 
of the DGLAP evolution with a very good accuracy; and (ii) the suppression of anti-DGLAP evolution of 
the transverse momenta of partons in the non- linear equation (see section 2.1 of this paper). 

It should be stressed that we reconstruct the form of the non-linear term in the modified BK equation 
assuming the probability interpretation of the parton cascade which was proven for the LO BK but has so 
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clear physical meaning that we believe that it will be correct in the NLO. However, the NLO corrections 
also lead to the transition of one dipole decaying into three dipoles, and induced by this transition term in 
the linear part of the equation. We neglected these terms since they induce the iV 3 term in the equation. 
This term is small in the vicinity of the critical line ( N oc as on it [31]), and only this kinematic region 
contribute to inclusive production. 

The failure to obtain the large corrections from the energy conservation suggests that the non- 
perturbative contribution is the reason for the small value of the NMF( nuclear modification factor). 
We wish to stress that the dipole rescatterings has been taken into account in Eq. ( [3. 36 ). The separate 
question is can we trust the kt factorization formula of Eq. (3.36) in the case of nucleus- nucleus collisions? 
The rigorous answer is no, since there is no proof of the kt factorization for this case. However, it is difficult 
to believe that the violation of the kt factorization will be responsible for the value of NMF. Indeed, in 
Eq. (HH) and Eq. (|1|) the main contribution stems from the vicinity of the saturation scale (the largest 
one in nucleus-nucleus scattering) and in the case of two scale problems: high px and the saturation scale 
the kt factorization has been proven in Ref. [29]. 

In our opinion the most interesting result of this solution is the appearance of the new scale Qf r = 
Q s exp(— 1/(20:5)), which characterize the temperature of the produced hadron ( T as Qf r ) . The produc- 
tion of dipoles with the sizes larger than 2/Qt r are suppressed. Much more work is needed to answer the 
question whether this scale is the artifact of our simplification or a general feature of energy conservation. 

As far as b dependence is concerned, we found that the BK equation without 6-dependence cannot 
guarantee an accuracy of calculation better than ±25%. We consider an important observation that in 
our simplified equation the entire 6-dependence can be absorbed in the non-perturbative behavior of the 
saturation scale. We hope that this fact can lead to a general scheme for taking into account the impact 
parameter dependence in the framework of high density QCD, however to rewrite the non-linear equation 
that absorbs the new non-perturbative scale, still needs more work. 
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